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Abstract 

We present results of a study of the muon decay in orbit (DIO) con- 
tribution to the signal region of muon - electron conversion. Electrons 
from DIO are the dominant source of background for muon - electron 
conversion experiments because the endpoint of DIO electrons is the 
same as the energy of electrons from elastic muon - electron conversion. 

The probability of DIO contribution to the signal region was con- 
sidered for a tracker with Gaussian resolution function and with a re- 
alistic resolution function obtained in the application of pattern recog- 
nition and momentum reconstruction Kalman filter based procedure 
to GEANT simulated DIO events. It is found that the existence of 
non Gaussian tails in the realistic resolution function does not lead to 
a significant increase in DIO contribution to the signal region. 

The probability of DIO contribution to the calorimeter signal was 
studied in dependence on the resolution, assuming a Gaussian resolu- 
tion function of calorimeter. In this study the geometrical acceptance 
played an important role, suppressing DIO contribution of the inter- 
mediate range electrons from muon decay in orbit. 
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Theoretical electron spectrum from muon decay in 
orbit 



Electron spectrum near the endpoint 

Electrons from muon decay in orbit are the dominant source of background 
for muon - electron conversion experiments. In the free decay of a muon at 
rest to an electron and two neutrinos, the electron's energy is at most half the 
muon rest energy, but in the decay of a bound muon the energy approaches 
that of the conversion electron, ~ 105 MeV, when the two neutrinos carry 
away little energy. 

General formulas obtained in [Ij and [2i|, describing the electron spec- 
trum from muon decay in orbit (DIO), are complicated. 

The results of numerical calculation of the electron spectrum for alu- 
minum and magnesium are presented in Table [TJ 

Note that in Table [T]the electron spectrum for aluminum is normalized 
to 1 but the electron spectrum for magnesium is not normalized. 

Near the endpoint the electron spectrum can be presented in an analyt- 
ical form found by Shanker [3]: 

N{E) = 10-21 . ( A)2 . ( A)5 . ff) + ^( A) + F{^)) (1) 

where 5 = Emax — E, 5i = — E — E'^ /2MA,and E is the elec- 
tron energy. The coefficients D, E, F depending on nuclear charge Z, were 
calculated for a wide range of elements in [3]. 

Eq.([T|) was obtained by neglecting the variation of the weak- interaction 
matrix element with energy. 

The maximum possible energy of the electron is: 

Emax = E^-^, E^ = m^-il-{a- Zf/2) 

where Ma is a nuclear mass, is a muon mass, a is the fine structure 
constant. In particular, E^nax is 104.963 MeV for aluminum. 

In order to apply Eq.([T]) to aluminum the coefficients D, E, F should be 
found for this element. This can be done, for example, by fitting numerical 
results presented in [3] by a polynomial of the 4th power (see Appendix A). 

As a result, for aluminum these coefficients are found to be 0.3575, 
0.9483, 2.2706, respectively. Note that these coefficients correspond to the 
normalization of the electron spectrum by the muon life time in aluminum 
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'^''Al Watanabe [4J 


""*M5 Herzog [2] 


Energy (MeV) 


N(E) (MeV-^) 




Energy (MeV) 


N(E) (MeV""^) 


5.0 


1.17-10-^ 




5.11 


1.249 • 10"^ 


10.0 


4.1 • 10"^ 




10.22 


4.382 • 10"^ 


15.0 


8.23 • 10"^ 




15.33 


8.814 • 10"^ 


20.0 


1.31 • 10"^ 




20.44 


1.41 


25.0 


1.85 • 10"^ 




25.55 


1.981 


30.0 


2.38 • 10"^ 




30.66 


2.549 


35.0 


2.87-10"^ 




35.77 


3.068 


40.0 


3.28 • 10"^ 




40.88 


3.487 


45.0 


3.52 • 10-^ 




45.99 


3.706 


50.0 


2.91 • 10"^ 




51.1 


2.683 


55.0 


3.67-10"^ 




56.21 


1.751 • 10"^ 


60.0 


1.41 • lO""* 




61.32 


6.056 • 10"^ 


65.0 


9.97-10-^ 




66.43 


4.44 • 10-4 


70.0 


1.11 • 10"*^ 




71.54 


5.008 • 10-^ 


75.0 


1.54 • lO" *" 




76.65 


6.846 • 10"^ 


80.0 


2.28 • 10-^ 




81.76 


9.704 • 10" '' 


85.0 


3.18 • 10"'^ 




86.87 


1.237 • 10" '' 


90.0 


3.54 • IQ-'^^ 




91.98 


1.158 • 10"^ 


95.0 


2.33 • 10"^^ 




97.09 


5.112 • lO"^'^ 


100.0 


3.58 • 10"^^ 




102.2 


1.632 • 10"^^ 



Table 1: Decay electron spectra (DIO) for '^'^Al and '^'^^Mg 



[5] tai = 0.864 lisec. The relation between the electron spectrum normal- 
ized by the muon life time in aluminum and the spectrum normalized to free 
muon decay is given by 

NiE) = NfreeiE)^ = Nfree{E)^^ « 0.4 • N free[E) 

where r/^ee = 2.2 jjLsec. Below N(E) will always refer to the spectrum 
normalized to free muon decay. 

According to Eq.( [T]) at E = 100 MeV the electron spectrum for alu- 
minum is N(E) = 1.42 X 10"^^ MeV~^. This number should be compared 
with N{E) = 1.43 x 10"^^ MeV~^, which is the result of numerical calcula- 
tions of 0] given in Table [T]but renormalized in order to take into account 
a muon life time in Al. On can see that the difference in these numbers is 
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less than 1%. At 95 MeV this approximation underestimates the spectrum 
by approximately 30 %. Remind that Watanabe's spectrum in Table [T]was 
normalized to free muon decay and calculated for energies below 100 MeV. 

Above E = 100 MeV the precision of Shanker's formula, obtained in 
phase space approximation by neglecting the variation of the matrix element 
with energy. Therefore the description of the spectrum is improved as one 
approaches the endpoint, because near the endpoint the DIO process is 
defined by the available phase space. Note that the numerical calculations 
[1] properly take into account relativistic electron wave functions and the 
effect of finite nuclear size on the wave functions. Thus we conclude that 
above 100 MeV Shanker's approximation ([T]) gives a good description of the 
electron spectrum with the precision better than 1%. 

It is important to emphasize that neglecting of the coefficients E and F 
at about 100 MeV leads to a significant underestimation (about 40%) of the 
electron spectrum N{E) = 1.0 x 10"^^ MeV~'^. 

In order to simplify analytical and numerical calculations and analysis 
of the results we decided to use, instead of Shanker's formula, a simple two 
terms approximation to describe the electron spectrum near the endpoint: 

N{E) = Co ■ {Emax — E)^ + Ci ■ {Emax — E)^ (2) 

For aluminum the coefficients Cq and Ci are determined by the following 
conditions: 1) It is required that at E = 100 MeV Eq.( [21) should give 
the same results N{E) = 1.43 • 10^^'^ MeV~^ as Watanabe's spectrum in 
Table [H renormalized in order to take into account the muon life time 
in Al. 2) Shanker's formula Eq.( [1]) is a good approximation above 100 
MeV, therefore at specific energy (we choose this energy as 104 MeV) near 
the endpoint the approximate formula ( [2|) should give the same result as 
Eq.( [5). These conditions give Co = 0.3699 x 10~^^ • MeV~^ and Ci = 
0.02132 X 10"^^ • MeV-'^ for aluminum. It was found that in this case in 
the range 100 - 104.5 MeV the spectra differ by less than 0.7%. 

At 100 MeV the difference between Shanker's approximation ( [T|) and 
two terms approximation ([2]) (and Watanabe's result) is about 0.7 %. At 95 
MeV these approximations underestimate the spectrum by approximately 
30 %. If one keeps only the leading term of Eq.( [2]) below 103.5 MeV the 
difference with Shanker's approximation is about 3 - 8 % 

Expected number of DIO events 

The expected number N of primary DIO events produced in aluminum target 
during the time of the experiment can be calculated as 



4 



^ — Ip ' ^^/p ' ^gate ' ^acc ' T ■ P (3) 

In this equation P is the probabihty of DIO contribution to a signal 
region above the threshold energy. By assuming that a proton flux Ip = 
4 X 10^^ /sec, muons stopping efficiency e^/p in the Al target per primary 
proton is 0.25%, Eacc = 0.2 is the setup overall acceptance, the time of 
experiment T = 10^ sec, and that the measured efficiency Sgate during the 
650 nsec window, extending from 700 nsec to 1350 nsec after the pulse, is 
50%: 

AT = 5 • 10^^ • P ■ Eacc- (4) 

Approximated theoretical electron spectrum 

Since the calorimeter energy cutoff can be as low as 80 MeV and the contri- 
bution from the intermediate part of the electron spectrum can be impor- 
tant. 

Using the numerical results from Table [Ufor aluminum, the theoretical 
electron spectrum can be approximated by a fit. 

Figure [1] shows the differential electron spectrum for muon decay in 
orbit in Al in linear scale. Points are the results of numerical calculations 
[4]. The solid line is an 8 parameter fit of the spectrum below 55 MeV in 
the form e-^^^^ where f(E) is a polynomial of power 7. The parameters of 
the fit are presented in Appendix B. 

Figure [2]shows the differential electron spectrum for muon decay in orbit 
in Al,a linear scale. The points are the results of numerical calculations [4j. 
The solid line is an 8 parameter fit of the spectrum in the range 55 - 100 MeV 
in the form e^^^^^ where g(E) is a polynomial of power 7. The parameters of 
the fit are presented in Appendix B. 

The region above 100 MeV is described by the two term approximation 
Eq.([2]) but its contribution in the case of low threshold energies is negligible. 

The probability to produce an electron of energy above Emin in muon 
decay in orbit in Al on the basis of the approximate theoretical electron 
spectrum is given by 

Emax 

P= j N{E)dE. (5) 

In this expression as it was mentioned above E^ax = 104.963 MeV , E 
is the true energy of electron (at this energy an electron is created in muon 
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Muon Decay in Orbit on AI27 
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Figure 1: Differential electron spectrum for muon decay in orbit in Al in 
a linear scale. Points are the results of numerical calculations. Solid line 
is a fit of the spectrum below 55 MeV. This spectrum is normalized to free 
muon decay. 

decay in orbit). 

Figure [3] shows the probability to produce an electron of energy above 
Emin in Al versus 

In the range 55 - 100 MeV the approximation discussed above was used, 
above 100 MeV the electron spectrum was taken in the form ([2]). 

Using the probability, which is given by Eq. ( [5]) , one can calculate the 
number of primary DIO events during the time of the experiment. It is 
important to emphasize that in this section we calculate the total number 
of DIO electrons produced in aluminum target and we do not need to include 
in this calculation a resolution function and overall acceptance. 

The number of DIO events, which is calculated using Eq.(l3|) with Sacc = 
1, is sharply increasing with decreasing Emin- In particular , for Emin = 100 
MeV - iV = 5.7 X 10^ for E^in = 95 MeV - TV = 6.2 x 10^ and for Emin = 
90 MeV - N = 1.39 x 10^. 

Table [2] shows the probability P and the number of DIO events N with 
an electron energy in the given energy range. 
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Muon Decay in Orbit on AI27 
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Figure 2: Differential electron spectrum for muon decay in orbit in Al in 
logarithmic scale. Points are the results of numerical calculations. Solid 
line is a fit of the spectrum in the range 55 - 100 MeV. This spectrum is 
normalized to free muon decay. 

Also in Table [21 the coefficients A and B are given, which are necessary 
for Monte Carlo simulation of DIO events. For this simulation it is assumed 
that each 5 MeV interval from Ei to E2 of the electron spectrum in the 
range 55 - 100 MeV can be approximated by 

N{E) = A ■ e-(-^-^i)/^ (6) 

It is important to note that in order to simulate a process with the 
threshold in measured energy E^ one has to take into account DIO elec- 
trons produced in target starting approximately at E^^ = E^ — 2a. In par- 
ticular, for Ejll = SOMeV and a = 5MeV it would correspond to E*^= 70 
MeV. According to Figure [3] and Eq. ([3]) in this case the number of primary 
DIO events would reach = 5.65 x 10^^ making event by event simula- 
tion unfeasible. We assume that present realistic threshold for calorimeter 
measured energy is about 90 MeV and the minimal true energy of DIO 
electrons is about 80 MeV. This corresponds to the number of DIO events 

= 1.13 X lO^*' which is about 50 times less than in the case of 80 MeV 
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DIO Integral Probability 



S 10"^ 

lio' 

Co 
S3 

10"^ 

10" 
10" 
10" 
10" 
10" 
10" 
10" 
10" 
in-i8 






















































































































































































































































..I. .1.1.. I.. 




,..,1.,.,...,.. 







70 75 80 85 90 95 100 105 110 
Minimal Energy, (MeV) 



Figure 3: Probability to produce an electron of energy above Emin in muon 
decay in orbit in Al 



Range 


80 - 85 


85 - 90 


90 - 95 


95 - 100 


> 100 


Prob. 


1.97 X 10'* 


2.62 X 10"^ 


2.64 X 10'^" 


1.23 X 10'" 


1.15 X 10'" 


N 


9.85 X 10^ 


1.31 X lO'* 


1.32 X 10* 


6.15 X 10^ 


5.7 X lO'' 


A 


9.12 X 10"^ 


1.272 X 10"^ 


1.416 X 10'"' 


9.32 X 10"^^ 




B 


2.53823 


2.27755 


1.83767 


1.19741 





Table 2: Expected number N of primary electrons produced in muon decay 
in orbit for different energy ranges (in MeV) . Coefficients A and B of Eq. Q 
are given. 



threshold in measured energy. 
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DIO events in a tracker 



Probability of DIO contribution 

The probability to have a signal from muon decay in orbit above a threshold 
is given by 

P= J cLEm j N{E)- f{EM-E)dE (7) 

-Emai + A 

where Em is the measured electron energy, E is the true energy at which 
an electron was emitted, A is the threshold energy measured from the end- 
point, N(E) is the electron spectrum of muon decay in orbit, and f is the 
resolution function of the tracker. 

For an ideal tracker the resolution function is the delta function and in 
this limit the probability P is given by 

Em ax 

P= f N{EM)dEM 

Em ax~\~^ 

For the high energy part of the electron spectrum, N{E) can be approx- 
imated by Eq.([2]), leading to the probability 

P = Co— -Ci— (8) 

in the case of the ideal detector. For the threshold energy A = 
-0.7 MeV, according to Eq.®, the probability P 6.9 • 10"^^. This 
equation is approximately valid for the range —5 MeV < A < 0. The lower 
limit is defined by the validity of the two term approximation ( [2|) for the 
electron spectrum. For positive A the probability P is for the ideal tracker. 

To deal with a realistic resolution function it is convenient to introduce 
a new variable y = Em — Em. Eq.(I7|). The region of integration in y,EM 
variables is given in Figure [H 

By interchanging the order of integration in Em^V and substituting 
Em = y + E the probability to have a signal from muon decay in orbit 
above a threshold takes the form 

P= j f{y)dy j N{E)dE (9) 

A Emax+A-y 
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Emax + A Em 
Figure 4: Domain of integration in y, Em variables. 



A second integral over the region in y above Emax + A was neglected 
because the resolution function at such y is expected to be extremely small. 

It follows immediately from this representation that the probability P is 
defined by the resolution function above the threshold energy A. 



Gaussian detector response function 

Let's assume that the detector response function f(y) is of the Gaussian 
form: 

fciy) = exp{ - ^ ) (10) 

In Figure [5] the probability P of muon decay in orbit as a function of 
threshold energy A is shown for different resolutions a. 

The two term approximation for the electron spectrum Eq.dJ]) was used 
for this plot. For small threshold energies A ~ the probability is sensi- 
tive to the resolution a and it is proportional to a^. At A = -1 MeV the 
probabilities for a = 0.4 MeV and a = 0.1 MeV differ by a factor 4. 

Note that Eq.Q defines the limiting probability for Figure [5] as the 
resolution a tends to 0. For the threshold energy A = —0.7 MeV and the 
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Figure 5: Probability of DIO contribution to a signal region as a function 
of the threshold energy A for different resolutions a. 



resolution a = 0.2 MeV the probability P « 2.1 • 10~^^ is about three times 
greater than the probability P ~ 6.9- 10^^^ in the case of the ideal detector. 

If the threshold energy A > — 5 MeV the probability P, given by Eq.([9|), 
can be calculated analytically by using the two term approximation Eq.([2]) 
for the electron spectrum. In this case 

The coefficients Iq , Ii are given by 



^-(15 + 45n^ + 15m* + u^)Erfc{^) - u{33 + Uu'^ + M*)e-"'/2^ (12) 



|u(105+105M2+21n^+M6)Sr/c(^)+(48+87n2+20n*+u^)e-"'/2 



(13) 
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where u = A/a, Erfc{z) = 1 — j= J dt 

In the case of interest, A < — cr, Eq. (jlip can be approximated by a form 
convenient for analysis 



0-6 ^1 

P = Coy (15 + 45n2 + 15n^ + u^) - Ci — n(105 + lOSu^ + 21n^ + u^) (14) 

The precision of this approximation is better than 1% if A < —a. In the 
hmit |A| >> a for negative A this equation reproduces Eq.dS]). 

Simulated detector response function 

The resolution function obtained as a result of simulation differs from the 
Gaussian form due to non Gaussian tails, which appear due to multiple scat- 
tering, radiation processes and a non-ideal reconstruction procedure. A pat- 
tern recognition and track reconstruction procedure based on the Kalman fil- 
ter technique developed in |6J was applied taking into account backgrounds, 
delta-rays and straw inefficiency. The distribution in the difference between 
the incident momentum {Pj^^) reconstructed by the forward Kalman filter 
procedure and the generated incident momentum {Pin) is shown in Figure [6] 
in linear (left) and logarithmic (right) scale. According to this distribution 
the intrinsic tracker resolution is <t = 0.189 MeV/c if one fits the distribution 
by a Gaussian in the range -0.3 to 0.7 MeV/c . It follows from Figure [6] 
that outside the range -0.4 to 0.4 MeV the distribution is non Gaussian. 

Using this histogram, normalized by the total area of the histogram, as 
the resolution function, the probability to have a signal from muon decay in 
orbit above a threshold can be presented in the form 



P = ^^ j N{E)dE (15) 

Emax+^—Pi 

where rij is the number in i*'' bin of the histogram, n = ni is the total 
number of events, pi = (p™" -|-p-""^)/2, p™" = A. 

In the two term approximation Eq.dJ]) for the electron spectrum this 
equation can be rewritten as 



n 



(16) 
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Figure 6: Distribution in the difference between the input reconstructed 
momentum and the generated input momentum with the background rate 
of 550 kHz, delta-ray and straw efficiency 97%. 

Figure [7] shows the probabihty of DIO contribution to the signal region 
as a function of the threshold energy A for a resolution function of the 
tracker obtained from the simulation (Figure [6|). 

One can see that in a logarithmic scale the probability decreases approx- 
imately linearly with an increase in threshold energy. This behavior differs 
from the sharp decrease of the probability in the case of Gaussian resolution 
function and is related to non Gaussian tails. The existence of non Gaussian 
tails leads to an increase in probability of DIO to contribute to the signal 
region but this additional contribution is not very significant. For example 
at the threshold energy of interest A = -0.7MeV P = 2.46 • 10~^^ in the 
case of simulated resolution function and P = 1.93- 10~^^ in the assumption 
of Gaussian resolution function with a = 0.189AIeV/c 

By using the two term approximation for the electron spectrum and the 
simulated resolution function, the expected number of primary DIO events 
in the tracker during the experiment with an Al target has been calculated 
and is presented in Table El 
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Figure 7: Probability of DIO contribution to a signal region as a function 
of the threshold energy A for a resolution function of the tracker obtained 
from the simulation. 



A, MeV 


-0.9 


-0.7 


-0.5 


-0.3 


-0.1 


N 


0.68 


0.25 


0.087 


0.03 


0.014 



Table 3: Expected number N of DIO events in a tracker in dependence on 
the threshold energy A. Realistic resolution function and overall acceptance 
are included. 



DIO events in the calorimeter 

A calorimeter geometrical acceptance (Figure [8]) plays an important role in 
suppression of low energy charged particles. 

One can see that due to a rapid fall-off the geometrical acceptance sup- 
presses the electron spectrum below 70 MeV, reducing significantly in this 
way the number of DIO events detected by the calorimeter above the thresh- 
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Figure 8: Geometrical acceptance of the calorimeter. 



old energy. 

Figure [9] shows a differential energy spectrum of DIO electrons multi- 
plied by the calorimeter geometrical acceptance. 

One can see that in the region below 70 MeV the rapid fall-off of the 
geometrical acceptance with decreasing electron energy almost compensates 
the rapid growth of DIO electron spectrum. In particular in the range from 
70 MeV to 60 MeV an increase in the electron spectrum by 2 orders of 
magnitude is reduced to a factor 2 due to the fall-off in the geometrical 
acceptance. 

The probability P to have a signal in the calorimeter from muon decay 
in orbit above a threshold is given by Eq.( [9]) with the electron spectrum 
multiplied by the geometrical acceptance Fgeom- 

+A 

P= I fG{y)dy I N{E) ■ Fg,om{E)dE (17) 

where, as above, y = Em — E, and Em is the measured electron energy, 
E is the true energy at which an electron was emitted. 

In Figure [10] the probability P of muon decay in orbit contributing to 
the calorimeter signal as a function of threshold energy A is shown. The 
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Figure 9: Differential energy spectrum of DIO electrons multiplied by the 
calorimeter geometrical acceptance. This spectrum is normalized to free 
muon decay. 



resolution function of the calorimeter is assumed to be Gaussian. 

It follows from this plot that the probability is not extremely sensitive 
to the resolution a, which is the consequence of fast decreasing geometrical 
acceptance below 70 MeV. For example, for the measured energy above 80 
MeV (the threshold energy A = — 25 MeV) even the change in a from 5 
MeV/c to 10 MeV/c increases the probability of DIO contribution only by 
a factor 5. According to this plot the probability P is a fast decreasing 
function of the threshold energy in the considered range because above 55 
MeV the electron spectrum is a steep function of energy. 

The expected number N of primary DIO events in calorimeter during 
the time of the experiment can be estimated by Eq.(j3]) with Eacc = 1- 

By using the fit for the electron spectrum described above and the Gaus- 
sian resolution function, the expected number of primary DIO events in 
calorimeter during the experiment with Al target is calculated and is pre- 
sented in Table [4] for a = 5 MeV/c and 8 MeV/c. 
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Figure 10: Probability of DIO contribution to the calorimeter signal as a 
function of the threshold energy A for different resolutions a. 



A, MeV 


-30 


-25 


-20 


-15 


-10 


N (a = 5 MeV) 


1.8 X 10"' 


5 X 10^ 


1.2 X 10^ 


2.5 X 10* 


4.3 X 10^ 


N (cr = 8 MeV) 


3.5 X 10^° 


1.4 X 10"' 


4.5 X 10^ 


1.3 X lO'' 


3 X 10* 



Table 4: Expected number N of DIO events in calorimeter in dependence 
on the threshold energy A. Convolution and geometrical acceptance are 
included. 



Muon conversion 

Muon conversion process is characterized by the appearance of a mono- 
energetic electron with energy E^°'^ = 104.963 MeV for aluminum. The 
rate of /ue - conversion R^^, is normalized by the muon capture rate: 
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where F^e is a width for muon conversion on nucleus and T^cap is a width 
for muon capture by nucleus. 
The total width is defined by 

^ fiTotal — ^fj,free ~l~ ^ficap 

or in terms of life times 

^— = ^ + ^ 

where for aluminum T^xotai = 0.864/isec and T^f^ee = 2.2/isec. 
The probability P^g for muon conversion is expressed as: 

P^ie = r^e/r^TotaZ = Rfie ' ^ ^icap ^iTotal = 0.6 • R^e 

The probability to have a signal from muon conversion on aluminum is 
given by: 

poo l-OO 

Pt^e' = P^^e ■ f{y)dy = 0.6 • R,,e ■ f{y)dy (18) 
The background to the signal ratio is expressed as: 

p-yp-^ = 0.103 • (10-iV^.e) • /r (y - A)^ f{y)dy/j^ fiy)dy 

where it was used that Co/(6 • P^e) = 0.103 • lO^^^/P^e- 

For a Gaussian resolution function this ratio can be approximated by 

Pow/Pt^' = 0-103 • (W-^VR^e) ■ • (15 + 45^.2 + is^x^ + ^6) 

where as above u = A/ a. The precision of this formula is better than 
2% for u < -2 . 

The expected number of registered muon conversion events can be cal- 
culated from Eqs.([5|).( fTB|) . For i?^e = 10~^^ sensitivity and setup overall 
acceptance = 20%, the expected number of muon conversion events is: 

iV = 4 • 10^3 X 2.5 • 10-3 X 10^ X 0.5 x 0.2 x 0.6 x lO'^^ = 6.0 
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Conclusion 



In this memo we studied the detection of electrons from muon decay in orbit. 
These electrons are the dominant source of background for muon - electron 
conversion experiments because the endpoint of DIO electrons is the same 
as the energy of electrons from elastic muon - electron conversion. 

It was found that near the endpoint (E ~ lOOMeV) the Shanker's formula 
[3] , obtained by neglecting the variation of the weak-interaction matrix ele- 
ments with energy in DIO process is in a good agreement with the results of 
numerical calculations [3j properly taking into account relativistic electron 
wave functions and the effect of finite nuclear size on the wave functions. 

It is important to note that in order to simulate a process with the 
threshold in measured energy one has to take into account DIO electrons 
produced in target starting approximately at E^^ = E^ — 2a. In particular, 
for E*J^ = 80MeF and a = 5MeV it would correspond to E^^= 70 MeV. 
In this case the number of primary DIO events would reach = 5.65 x 
10^^ making event by event simulation unfeasible. We assume that present 
realistic threshold for calorimeter measured energy is about 90 MeV and the 
minimal true energy of DIO electrons is about 80 MeV. This corresponds 
to the number of DIO events = 1.13 x 10^*^ which is about 50 times less 
than in the case of 80 MeV threshold in measured energy. 

The probability of DIO contribution to a signal region was considered 
for the tracker with Gaussian resolution function and with the realistic res- 
olution function obtained in the application of pattern recognition and mo- 
mentum reconstruction Kalman filter based procedure to GEANT simulated 
DIO events. 

It was found that non Gaussian tails in the simulated resolution func- 
tion do not lead to a significant increase in DIO contribution to the signal 
region. The expected number of detected DIO events during the time of the 
experiment for realistic resolution function was calculated to be about 0.25 
if the threshold energy A = —0.7 MeV and the overall acceptance is 20%. 

It was shown that the probability of DIO contribution is very sensitive 
to the threshold energy and proportional to in the limit |A| ^ o" where 
a is the tracker resolution. 

The probability of DIO contribution to the calorimeter signal was stud- 
ied in dependence on the resolution, assuming a Gaussian resolution func- 
tion of the calorimeter. In this study the geometrical acceptance played 
an important role, suppressing DIO contribution of the intermediate range 
electrons from muon decay in orbit. It was found that the probability of 
DIO contribution is not extremely sensitive to the resolution a, which is the 
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consequence of the fast decreasing geometrical acceptance below 70 MeV. 
For example, for the measured energy above 80 MeV (the threshold energy 
A > —25 MeV) even the change in a from 5 MeV/c to 10 MeV/c increases 
the probability of DIO contribution only by a factor 5. 

The expected number N of detected DIO events in calorimeter was esti- 
mated by using the approximation of the electron spectrum in the range 55 - 
100 MeV. For the measured energy above 80 MeV and assuming a Gaussian 
resolution function of the calorimeter with a = 5 MeV it was found that 
during the time of the experiment N = 5 x 10^. 

We wish to thank A. Mincer and P.Nemethy for fruitful discussions and 
helpful remarks. 

References 

[1] P.Hanggi, R.D.Viollier, U.Raff, and K.Alder, Phys.Lett. 51B, 119 
(1974) 

[2] F.Herzog and K.Alder, Helv.Phys.Acta 53, 53 (1980) 
[3] O.Shanker, Phys.Rev. D25, 1847 (1982) 

[4] R.Watanabe et al.. Atomic and Nuclear Data Tab., 54, 165 (1993) 
[5] T.Suzuki, D.F.Measday and J.P.Roalsvig, Phys.Rev. C35, 2212 (1987) 
[6] R.Djilkibaev and R.Konophch, |hep-ex/03T2022| (2003) 



20 



Appendix A 

Coefficients of Shanker's expansion 



The coefficient D, E, F oi Shanker's expansion for different materials can 
be determined by the fitting of numerical values of D, E, F presented in [3] 
by 4th power polynomial. 

Figure [TT] shows the results of the fit for D coefficient. 




Figure 11: Numerical values of D coefficient. Solid line is a result of fit by 
4th power polynomial / = pi + p2 * x + p3 * + p4 * + p5 * x^. 

In the same way E and F coefficients can be found. 
Figures [12] and [T3] show the results of the fit for D and E coefficients, 
respectively. 

Table [5] represents D, E, F coefficients for different materials. 



21 




Figure 12: Numerical values of E coefficient. Solid line is a result of fit by 
4tli power polynomial f = pi + p2 * x + p3 * x'^ + pA* + p5 * x^. 

Appendix B 

Fit of electron spectrum 

Below 55 MeV the electron spectrum Nf^eeiE) can be described by an 8 
parameter fit in the form e^^^^ where f(E) is a polynomial of power 7: 

f{E) = h-E'' + h-E'' + h-E'> + U-E^ + h-E^ + f2-E^ + h-E + h (19) 

For aluminum the fit, using numerical results for the electron spectrum 
[1], gives: 

/7 = -4.1321585902527747 x 10"^°; 
/e = 7.26144332670667 x 10"^ 
h = -5.0391208057229 x 10"^ 
/4 = 1.720029756203247 x lO""^; 
h = -2.797753809879399 x lO-^; 
/2 = 1.003520948296598 x 10^^. 
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Figure 13: Numerical values of F coefficient. Solid line is a result of fit by 
4th power polynomial / = pi + p2 * x + p3 * x'^ + pA* + p5 * x'^. 

fi = 0.35027124751754113; 
/o = -8.495230158021982; 

Above 55 MeV the electron spectrum Nfree{E) can be described by an 
8 parameter fit in the form e^^^^ where g(E) is a polynomial of power 7: 



For aluminum the fit, using numerical results for the electron spectrum 
[4], gives: 

gj = 5.882860554682616 x 10"^°; 
56 = -3.311494207831439 x lO"'^; 

05 = 7.90482092857202 x 10"^; 
54 = -1.0375709354278761 x lO^^; 
53 = 0.808726145419064; 
£,2 = -37.4179819613308; 
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Table 5: E, F coefficients for different materials 











z 


D 


E 


F 


10 


0.169990093 


0.281679779 


1.68660045 


11 


0.221416339 


0.487800866 


1.8169477 


12 


0.282654375 


0.704645038 


1.99828386 


13 


0.357468009 


0.948250234 


2.27059889 


14 


0.449176699 


1.2330687 


2.6695714 


15 


0.560655653 


1.57196677 


3.22657013 


16 


0.694335938 


1.97622538 


3.96865273 


17 


0.852204204 


2.45553946 


4.91856623 


18 


1.03580284 


3.01801825 


6.09474707 


19 


1.24623013 


3.67018533 


7.51132059 


20 


1.48413992 


4.41697836 


9.17810249 


21 


1.74974191 


5.26174974 


11.1005974 


22 


2.04280138 


6.20626593 


13.2799997 


23 


2.36263967 


7.25070715 


15.713191 


24 


2.70813346 


8.39366913 


18.3927441 



gi = 950.4758914390899; 
go = -10212.91359983421; 
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